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We study the dynamics of vortex-antivortex (VA) pairs in an infinitely thin ferromagnetic film 
with easy-plane anisotropy. These are localized excitations with finite energy that are characterized 
by a topological (skyrmion) number M = 0, ±1. Topologically trivial {N = 0) VA pairs undergo 
Kelvin motion analogous to that encountered in fluid dynamics. In contrast, topologically nontrivial 
(Af — ±1) VA pairs perform rotational motion around a fixed guiding center. We present the results 
of a detailed study in both cases and further demonstrate that in the presence of dissipation a 
rotating J\f = ±1 VA pair shrinks to a point and is annihilated, due to the discreteness of the 
lattice, thus leading to a "topologically forbidden" AJV = I process. We argue that the latter 
process underlies the experimentally observed vortex core switching whereby the polarity of a single 
vortex is reversed after collision with an Af = VA pair created by a burst of an applied alternating 
magnetic field. 



I. INTRODUCTION 



The best known examples of topological magnetic solitons are magnetic bubbles or skyrmions observed 
in abundance in ferromagnetic films with easy- axis anisotropy [H, 0j- The experimental situation is less 
clear in the case of ferromagnets with easy-plane anisotropy. The relevant topological structures are 
theoretically predicted to be half skyrmions or vortices, with a logarithmically divergent energy that 
may inhibit production of an isolated vortex on an infinite film. Thus the early studies of magnetic 
vortex dynamics had been mostly theoretical [H, 0] drawing on various analogies with related work on 
ferromagnetic bubbles Q , with vortex dynamics in classical fluids and superfluids [H, 0, 0] , as well as 
with the dynamics of interacting electric charges in a uniform magnetic field. 

The situation has changed dramatically in recent years. It has been realized that a disc-shaped mag- 
netic element, with a diameter of a few hundred of nanometers, provides an excellent geometry for the 
realization of a magnetic vortex configuration. In particular, the exchange energy is finite on a finite ele- 
ment while the magnetostatic field vanishes everywhere except at the vortex core. As a result, the vortex 
is actually the lowest energy magnetic state in a disc-shaped element. In other words, interest in the 
vortex stems from the fact that it is a nontrivial magnetic state which can, nevertheless, be spontaneously 
created in magnetic elements 

It is then natural to ask whether nontrivial magnetic states other than the single vortex may play an 
important role in the dynamics of magnetic elements 1^, 11 1 • An answer to this question comes from a 
somewhat unlikely direction. Recent experiments have shown a peculiar dynamical behavior of vortices 
and magnetic domain walls when these are probed by external magnetic fields. Vortices may switch their 
polarity under the influence of a very weak external magnetic field of the order of a few mT [T3 . fl^ . 
The same switching phenomenon was observed by passing an a.c. electrical current through a magnetic 
disc [14|. Since the polarity of the vortex contributes to its topological structure, the switching process 
clearly implies a discontinuous (topologically forbidden) change of the magnetic configuration. This is 
certainly a surprise especially because the external field is rather weak. The key to this phenomenon is 
the appearance of vortex pairs which are spontaneously created in the vicinity of existing vortices [l3L ITEj . 
The creation of topological excitations (vortex pairs) by alternating external fields had been anticipated 
by an early study based on collective coordinates [lq . 

In this paper we study vortex-antivortex pairs (VA pairs) which are nontrivial magnetic states that play 
an important role in the dynamics of magnetic elements. However, unlike a single vortex, a VA pair is a 
localized object whose energy remains finite even on an infinite film. It is then reasonable to expect that 
the essential features of the dynamics of VA pairs can be understood in the infinite-film approximation 
which is adopted in the following. A brief summary of the relevant dynamical equations and related 
topological structures is given in Section [TTJ In Section IIIII we study a VA pair in which the vortex and 
the antivortex carry the same polarity. Such a pair is shown to undergo translational Kelvin motion 
analogous to that observed in fluid dynamics [17j. In Section QV] we study a VA pair in which the vortex 
and the antivortex carry opposite polarities. Such a pair is shown to behave as a rotating vortex dipole 
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[l8| because its topological structure is substantially different than that of the pair in Kelvin motion. 

The preceding results are combined in Section [V] to demonstrate that a rotating vortex dipole may be 
annihilated by a quasi-continuous process in spite of its nontrivial topological structure. In particular, 
no energy barrier has to be overcome in contrast to the usual expectations for topological solitons. This 
opens the possibility for switching mechanisms between topologically distinct states in ferromagnets. 
The possibility to change the topological structure leads to a dramatic change in the magnetization 
dynamics as a VA pair is created or annihilated. Such a pair annihilation process lies in the heart of the 
counter- intuitive vortex polarity switching event that was observed in magnetic elements 12 , 13|, 3 ■ 



Some concluding remarks are summarized in Section IVII Finally, an Appendix is devoted to a brief 
description of the dynamics of interacting electric charges in a uniform magnetic field, which also exhibits 
most of the peculiar features of the dynamics of VA pairs. 

II. THE MODEL 

A ferromagnet is characterized by the magnetization m = (mi, m2, ma) measured in units of the 
constant saturation magnetization M s . Hence m is a vector field of unit length, m 2 — m 2 + m 2 . + m 2 = 1, 
but is otherwise a nontrivial function of position and time m — m(r,t) that satisfies the rationalized 
Landau-Lifshitz (LL) equation 

^ = m X /. f = Am-qm 3 e 3 , m 2 = 1. (1) 



Here distances are measured in units of the exchange length £ cx = yj A/2-kM 2 , where A is the exchange 
constant, and the unit of time is tq = 1/(Att^M s ) where 7 is the gyromagnetic ratio. Typical values are 
f cx ~ 5nm and To ~ lOps which set the scales for the phenomena described by Eq. ([1]). To complete the 
description of the LL equation we note that we consider ferromagnetic materials with uniaxial anisotropy. 
Then (=3 in Eq. {T]) is a unit vector along the symmetry axis and the dimensionless parameter q = K/2irM 2 , 
where K is an anisotropy constant, measures the strength of anisotropy. In particular, q is taken to be 
positive throughout this paper, a choice that corresponds to easy-plane ferromagnets. 

An important omission in Eq. |T]) is the demagnetizing field produced by the magnetization itself 0, Q • 
However, in the limit of a very thin film, the effect of the de mag netizing field is thought to amount to 
a simple additive renormalization of the anisotropy constant [191 ]. Also note that we may perform the 
rescalings ^fqr — ► r and qt — > t which further renormalize the units of space and time discussed earlier 
and lead to a completely rationalized LL equation where we may set q = 1 without loss of generality. 
With this understanding, all calculations presented in this paper are based on a two-dimensional (2D) 
restriction of Eq. (JTJ) , i.e., r = (x,y) and A = d 2 /dx 2 + d 2 /dy 2 , while q is set equal to unity without 
further notice. 

The effective field / in Eq. ([1]) may be derived from a variational argument: 

/ = -||; E = l - J [(Vm) 2 + m%] dxdy (2) 

where E is the conserved energy functional. A standard Hamiltonian form is obtained by resolving the 
constraint m 2 = 1 through, say, the spherical parameterization 

m 1 + i TO2 = sin e 4 *, m 3 = cos O. (3) 

The LL equation is then written as 

d<$> _SE OH _ 8E 

In^m' ~at ~ ( ' 

where n = cos0 is the canonical momentum conjugate to the azimuthal angle $. Taking into account 
the specific form of the energy in Eq. j2]) or 

E =\j [(V©) 2 + sin2 © (V$) 2 + cos 2 6] dxdy (5) 
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the Hamilton equations Q yield 

sinG^- = A9 + [1 - (V$) 2 ] cosG sin 6, 
BP) 

sin6— = -V ■ (sin 2 e V*). (6) 
Another useful parameterization is obtained through the stereographic variable f2: 

mi+j?7i2 . 2f2 1 — Ofi 
it = ; mi + i = =— , m<? = =— (7) 

l + ma l + asi l + nn w 

in terms of which the LL equation reads 

+ An + 1^2^ n = — 5L- (vn • vn) (8) 
9t 1 + 1 + nn 

where f2 is the complex conjugate of il. Equations ([I]),©, and ijHJ) are three equivalent versions of the 
LL equation and may be used at convenience depending on the specific calculation considered. 

A key quantity for describing both topological and dynamical properties of the 2D LL equation is the 
local topological vorticity 7 = j(x,y,t) defined from [20L Ell] : 

7 = e a /3 d a Ildf3$ = e af3 sin Q dpOd a $ = ^e a0 (d m x d a m) ■ m, (9) 

where the usual summation convention is invoked for the repeated indices a and /3, which take over two 
distinct values corresponding to the two spatial coordinates x and y, and e a p is the 2D antisymmetric 
tensor. In particular, one may consider the total topological vorticity Y and the Pontryagin index or 
skyrmion number J\f defined from 

r = j jdxdy, Af = (10) 

A naive partial integration using Eq. ([9]) yields T = = N for all magnetic configurations for which such 
an integration is permissible. However, nonvanishing values for L and Af are possible and are topologically 
quantized. Specifically, for field configurations that approach a constant (uniform) magnetization at 
spatial infinity, the skyrmion number Af is quantized according to Af = 0, ±1, ±2, . . .. Half integer values 
are also possible in the case of field configurations with more complicated structure at infinity such as 
half skyrmions or vortices (see below). 

The local topological vorticity 7 is also important for an unambiguous definition of conservation laws 
in the LL equation. Hence the linear momentum (impulse) P = (P x ,P y ) is defined from 



P, = 



Jyydxdy, P y = J xjdxdy, (11) 



while the angular momentum (impulse) is given by 



p z jdxdy (12) 



where p 2 = x 2 + y 2 ■ Since detailed discussions of these conservation laws have already appeared in the 
literature 17, 2^, 21 1 we simply note here that analogous conservation laws were defined as moments of 



ordinary vorticity in fluid dynamics 0, 0] • 

We first search for static (time independent) solutions of the LL equation which may be obtained by 
omitting time derivatives in Eq. ([6|) and further introducing the axially symmetric ansatz Q = 9(p) and 
$ = k(4> — 4>o), where p and <j> are the usual cylindrical coordinates [x = pcoscj), y — psin0), k = ±1 
will be referred to as the vortex number, and 4>o is an arbitrary constant phase reflecting the azimuthal 
invariance. The resulting ordinary differential equation for the amplitude 6 — 8(p) is solved numerically 
with standard boundary condition 6{p — * 00) = 7r/2 and the result is shown in Fig.[T] The corresponding 
magnetization is then given by 

mi +im 2 = smOe^-^, m 3 = A cos (9, (13) 
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FIG. 1: Profile of a single vortex calculated numerically. The complete solution is given by Eq. (flU 



where A = ±1 will be called the polarity. The total energy is accordingly reduced to 

1 



E = 



06 

dp 



cos z e 



(2irpdp). 



(14) 



where the centrifugal (second) term is logarithmically divergent for the assumed boundary condition at 
spatial infinity. However, the anisotropy energy given by the last term is finite and is actually predicted 
to be 



Ea = \j C0S2 9 = \ 



(15) 



by a careful derivation of a suitable virial relation [221 ]. 
number M are calculated from Eq. (fit))) to be 

T = -2ttkA , TV 



Finally, the total vorticity T and the skyrmion 



(16) 



where the vortex number k = ±1 and the polarity A = ±1 may be taken in any combination. Thus, we 
must consider four possibilities; namely, a vortex that comes in two varieties (k = 1, A = ±1) and thus 
J\f = =Fl/2, and an antivortex which also comes in two varieties (k — —1, A = ±1) and thus Af = ±1/2. 
In all cases the calculated static solution is a topological soliton with half integer skyrmion number, in 
contrast to ordinary skyrmions (such as magnetic bubbles) which carry integer Af. 

Our main aim in the continuation of this paper is to search for nontrivial solutions that may combine 
a vortex and an antivortex (VA pair) in a way that the total energy is finite. We leave aside for the 
moment the LL equation and construct model VA pairs in terms of the basic single-vortex configuration 
described above. This is easily accomplished by invoking the stereographic variable of Eq. to write 
a single (k, A) vortex located at the origin of coordinates as 



n = 



shi9 



iK,(<j>-<j>o) 



1 + A cos ( 



(17) 



where 9 = 9{p) is the vortex profile taken from Fig. [T]or simply the model profile defined from cos, 9 = 
1 / cosh p and sin 9 = tanh p. We now produce two replicas of the basic vortex to describe a pair of vortices 
by the product ansatz 



(18) 



where Q\ is configuration (|17p applied for (k, A) = (ki,Ai) and the origin displaced to, say, (x,y) — 
(— d/2,0) while f2 2 is a (k 2 ,A 2 ) vortex located around (x,y) = (d/2,0). The skyrmion number of this 
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configuration is given by 



Af = -- (k x Ai + k 2 A 2 ) 



(19) 



and its energy is finite only if we restrict attention to vortex-antivortex (VA) pairs (k\ = — fv 2 ). For 
dcfiniteness we choose k± = 1 and k 2 = — 1 and thus the skyrmion number 



depends on the polarities Ai and A 2 . Now a VA pair with equal polarities (Ai = A 2 = ±1) is topologically 
trivial (Af — 0). In contrast, a VA pair with opposite polarities is topologically equivalent to a skyrmion 
(Af = 1, for Ai = —1 and A 2 = 1) or an antiskyrmion (Af = —1, for Ai = 1 and A 2 = —1). In all three cases 
configuration (|18p carries finite energy because its overall phase cancels out at spatial infinity where the 
magnetization approaches the uniform configuration m — (1,0,0) modulo an overall azimuthal rotation 
which depends on the choice of individual phases 0o in the ansatz (|18[) . This explains, in particular, why 
VA pairs are characterized by an integer skyrmion number (Af — 0,±1). 

Needless to say, the VA pairs constructed above are not solutions of the LL equation. Yet an interesting 
picture arises when configuration (|18|) is used as an initial condition in the complete equation ((8|). A 
topologically trivial (Af — 0) VA pair undergoes Kelvin motion in which the vortex and the antivortex 
initially located at a relative distance d along the x axis move in parallel along the y axis with nearly 
constant velocity. In contrast, a topologically nontrivial (Af = ±1) VA pair undergoes rotational motion 
around a fixed guiding center at nearly constant angular velocity. In both cases the main trajectories are 
decorated by Larmor-type oscillations [23| which are tamed when the relative distance between the vortex 
and the antivortex is large. The effect of the polarity of vortex pairs has been studied experimentally in 
patterned ferromagnetic ellipses [24^ . Two vortices were created an ellipse and it was found that their 
dynamics depended on their relative polarity, as is indicated by Eq. (j2"0|) . 

In the following two sections we examine the two cases in turn. In particular, we aim at constructing 
true steady-state solitary waves that describe VA pairs in pure Kelvin motion for V = and pure 
rotational motion for Af = ±1. 



As the title of this section suggests, VA pairs in Kelvin motion were originally studied in the context 
of ordinary fluid dynamics @, 0j. A further analogy exists with the 2D motion of an electron-positron 
pair interacting via the Coulomb potential and placed in a uniform magnetic field perpendicular to the 
plane. If the electron and the positron are initially at rest their guiding centers will move along two 
parallel straight lines while the actual trajectories will display the familiar Larmor oscillations. However, 
when the electron and the positron are given a common initial velocity such that the Coulomb force is 
exactly balanced by the magnetic force, both the guiding centers and the actual positions of the charges 
will move steadily along parallel lines, even though the two sets of trajectories do not coincide. The 
resulting special configuration may be thought of as a peculiar electron-positron bound state in steady 
translational motion (see our Appendix). 

It is thus reasonable to expect that a solitary wave exists in a 2D easy-plane ferromagnet which 
describes a VA pair that proceeds rigidly (without Larmor oscillations) in a direction perpendicular to 
the line connecting the vortex and the antivortex, probably because the mutual force is exactly balanced 
by a topological "Magnus force" . The actual existence of such a solitary wave was established in Ref. [13] 
whose main result is briefly reviewed in the remainder of this section. 

As it turns out, when the relative distance is large, the sought after solitary wave resembles in its 
gross features the model VA pair of Eq. (fT8|) applied for, say, Ki = —k 2 = 1 and Ai = A 2 = 1 (thus 
Af = 0). We may then invoke this model to motivate some important asymptotic results valid for large 
d. For instance, the local topological vorticity 7 is then peaked around the two points (x, y) = (—d/2, 0) 
and (x, y) = (d/2,0) with weights — 2ir and 27r respectively, corresponding to the total vorticities T of 
the individual vortex and antivortex. Then the impulse defined from Eq. (|11[) yields P x — 0, thanks to 
reflexion symmetry, while P = P y ~ 2nd. One may also invoke the Derrick-like scaling relation applied 
to the extended energy functional F = E — vP: 



Af = --(Ai-A 2 ) 



(20) 



III. KELVIN MOTION 




(21) 
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FIG. 2: Snapshot of a topologically trivial (A/ - = 0) VA pair in Kelvin motion, illustrated through the (mi, 7712) 
projection of the magnetization m = (mi, 7712, 7713) on the plane of the film (left panel) as well as level contours 
of m3 (right panel). The contour levels 7713 = 0.1,0.3,0.5,0.7,0.9 are shown. The outer curve corresponds to 
7773 = 0.1, while the two smallest circles are the contours for 777,3 = 0.9 and surround the two vortex centers. The 
pair moves along the y axis with velocity v = 0.2 at a calculated relative distance d — 4.97 (along the x axis), 
energy E = 20, and impulse P = 30. 



where E a is the total anisotropy energy of the VA pair. For large d this energy approaches the sum of 
the anisotropy energies of the individual vortex and antivortex, each given by 7r/2 in view of the virial 
relation (fl5|) . Hence, J 777,3 dxdy ~ 2(ir/2 + n/2) = 2ir and vP ~ 2tt. To summarize, 

P~2ird, vP~2ir, v . (22) 

a 



One may further consider the familiar group-velocity relation 

_ dE 
V ~dP 



(23) 



in which we insert the estimate v ~ 2n/P to obtain an elementary differential equation for E whose 
integral is 

E w 2tt ln(P/P ) (24) 

where Pq is an integration constant that cannot be fixed by the present leading-order argument. Never- 
theless, Eq. (|24[) provides the essence of the energy-momentum dispersion at large relative distance d or 
small velocity v ~ l/d and hence large momentum P ~ 2?td. 

To obtain an accurate numerical solution it is convenient to work with the stereographic variable f2 of 
Eq. ([7]). Then the LL Eq. ([8]) restricted to a solitary wave in rigid motion with constant velocity v along, 
say, the y axis reads 

- iv ^l + An + LzIB- n = — 5L- (vn • vn) (25) 
d y 1 + nn 1 + nn y ' ' 

and is supplemented by the boundary condition Q — > 1 at spatial infinity where the magnetization 
approaches a uniform configuration. Once a solution = £l(x, y\ v) of Eq. f|25|) is obtained for a specific 
value of the velocity v, the sought after solitary wave is given by f2(x, y — vt; v). 

Equation (|25p was solved numerically via a Newton-Raphson iterative algorithm for velocities in the 
range 0.1 < v < 0.99. Note that v = 1 is the familiar magnon velocity (in rationalized units) and provides 
an upper bound for the existence of a solitary wave in rigid motion. On the contrary, there is no lower 
bound for the velocity - the restriction v > 0.1 was dictated only by numerical expedience. Now, the 
calculated solitary wave is illustrated in Fig. [2] for the relatively low velocity v — 0.2 and does indeed 
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FIG. 3: Solitary wave in Kelvin motion with v = 0.95 along the y axis (conventions as in Fig. Note that the 
vortex-antivortex character is lost (d = 0) and the wave is a lump with no apparent topological features. The 
calculated energy and impulse are E = 18 and P — 17. 




W 20 




FIG. 4: Energy E and impulse P as functions of velocity v (left panel) and E vs P dispersion (right panel) for a 
solitary wave in Kelvin motion. The dotted line was calculated from the asymptotic dispersion (|26[1 applied for 
Qo = 4.7. 



describe a VA pair (with J\f = Q) at a relative distance d = 4.97 ~ 1/v, as anticipated by the discussion of 
model VA pairs in Section HI1 and the heuristic asymptotic analysis earlier in this section. But we are now 
in a position to carry out an accurate calculation practically throughout the allowed range < v < 1, to 
discover that there exists a characteristic velocity vq ~ 0.78 above which the vortex-antivortex character 
is lost (d ~ 0) and the solitary wave becomes a lump with no apparent topological features, as illustrated 
in Fig. for v = 0.95. 

The existence of a characteristic velocity v becomes apparent when we calculate energy E and impulse 
P as functions of v, as shown in Fig. [5] Note that both E and P develop a minimum at a common velocity 
v = vq = 0.78. As a result, the energy vs impulse dispersion shown in Fig. [4] develops a cusp at a point 
(Eq, Pq) that corresponds to the values of E and P at v = vq. Thus the calculated family of solitary waves 
consists of two branches. Branch I consists of VA pairs that propagate with velocities in the range v < vq. 
The corresponding branch in the dispersion of Fig. |4] approaches the asymptotic dispersion ()24|) for large 
P (or v — * 0). Indeed, an excellent fit of the data is obtained by Eq. if we choose the subleading 
constant according to InPo ~ 1/4- Branch II consists of lumps with no apparent topological features 
which propagate with velocities in the range vq < v < 1. The corresponding branch in the dispersion of 
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Fig. [4] is accurately described in the asymptotic (v — > 1) region, where P again becomes large, by 

E = p ( 1 + $i + ~) (28) 

with Q w 4.7. Actually this asymptotic dispersion can be derived by showing that in the limit v — > 1 
the 2D Landau-Lifshitz equation reduces to what is a modified Kadomtsev-Petviashvili (KP) equation 

E2|. 

A cusp in the energy vs impulse dispersion occurred previously in a calculation of vortex rings in a 
model superfluid by Jones and Roberts [251 ]. The same authors together with Putterman later argued 
that the solitary waves that correspond to branch II are actually unstable [26| . Interestingly, a dispersion 
with a cusp appears also in the much simpler problem of electron-positron motion in a uniform magnetic 
field, where the motion that corresponds to branch II is also unstable (see our Appendix). Hence, while 
a stability analysis has not yet been carried out for the solitary waves described in this section, it is 
reasonable to expect that the lumps of branch II may be unstable. But there is every reason to believe 
that the Kelvin motion of the VA pairs of branch I is indeed stable. 



IV. ROTATIONAL MOTION 



The possibility of topologically nontrivial (M — ±1) VA pairs in steady rotational motion was recently 
examined by one of us [181 ] . Again, one may invoke the model VA pair of Section [ill to understand 
some important features of the rotational motion at large distance d. For definiteness we consider a VA 
pair defined by Eq. (fTg|) with K\ = —k 2 = 1 and Ai = — A 2 = —1, thus J\f = 1. The local topological 
vorticity 7 is then peaked around the positions of the vortex and the antivortex, now with weight equal 
to 2ir in both cases. Then, for large d, the angular momentum defined by Eq. (fT2"j) is estimated to be 



L <~ i2(27r)(|) 2 = jd 2 . One may also consider the Derrick- like scaling relation applied to the extended 
energy functional F = E — ujL: 

ujL = — J m 2 dxdy = E a (27) 

where E a is the total anisotropy energy of a VA pair which approaches asymptotically B a ~(| + 5)=7r, 
and hence loL — ir, because the anisotropy energy of a single vortex is equal to ir/2 according to Eq. (p~5|) . 
To summarize, 

L~|d 2 , uiL-n, (28) 



One may further employ the familiar relation 



dE , . 

U = dL < 29 > 



in which we insert the estimate uj ~ tt/L to obtain an elementary differential equation for E whose 
integral is 

E «7rm(£/L ) (30) 

where Lo is an integration constant that cannot be fixed by the present leading-order argument. Nev- 
ertheless, Eq. (|30p provides the essence of the energy vs angular momentum dispersion for large relative 



distance d or small angular frequency u> ~ 2 / d 2 and hence large angular momentum L ~ ? d 2 . 

While the preceding heuristic asymptotic analysis is very useful for understanding some basic aspects of 
a rotating VA pair, it does not give us any clue concerning the fate of the pair at small vortex-antivortex 
separation. In principle, such a question could be settled by solving numerically the analog of Eq. (|25p 
for a solitary wave rotating at constant angular frequency lo: 

1 - Tin 20 

iuj e a g x a dgfl + Afl H = — (VO-Vfi). (31) 

A numerical solution could again be attempted by an iterative Newton-Raphson algorithm. Actually, in 
this case, we found it more convenient to employ a relaxation algorithm to derive approximate numerical 
solutions as stationary points of the extended energy functional F = E — ujL [lfl • 



FIG. 5: Snapshot of a topologically nontrivial (TV = 1) VA pair in rotational motion (conventions as in Fig. [2] 
with solid lines in the right panel corresponding to positive values of 7713 and dashed lines to negative ones). The 
pair rotates around a fixed guiding center taken at the origin of coordinates, with angular velocity u) — 0.06 and 
calculated relative distance d — 5.3, energy E = 21, and angular impulse L = 64. 
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FIG. 6: J\f = 1 VA pair in rotational motion with angular velocity u) — 0.18 and calculated relative distance 
d = 1.7, energy E — 15, and angular impulse L = 11. The only difference from Fig. [S] is that the overall size of 
the pair is now reduced. 



The calculated configuration for uj = 0.06 is shown in Fig. [5] and does indeed correspond to a topolog- 
ically nontrivial (AT — 1) rotating VA pair consisting of a vortex with negative polarity (k, A) = (1, —1) 
and an antivortex with positive polarity (re, A) = (—1,1), as anticipated by the general discussion of 
Section [TTJ In particular, for this relatively small value of w, the calculated distance d in the rotating 
pair is relatively large (d = 5.3), so is the angular momentum [L = 64); in rough agreement with the 
asymptotic estimates of Eq. p8]l . 

The calculation was repeated for a larger value of angular velocity (ui — 0.18) to find that both the 
relative distance [d = 1.7) and the angular momentum [L — 11) are reduced to smaller values. But the 
general structure of the solution shown in Fig. [6] for lo = 0.18 remains basically the same as that for 
u) = 0.06, except that the overall size of the VA pair is reduced. For yet larger values of u> the angular 
momentum L tends to vanish. This trend is apparent in Fig. [7] which illustrates the dependence of E 
and angular momentum L as functions of angular velocity u), as well as the E vs L dispersion. For large 
values of L the above dispersion exhibits logarithmic dependence, as predicted by the asymptotic result 
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FIG. 7: Energy E and angular impulse L as functions of angular velocity u> (left panel) and E vs L dispersion 
(right panel) for an N = 1 VA pair in rotational motion. 



of Eq. (|30j) . But the most important feature of the calculated dispersion is that energy approaches the 
finite value E = Att as L — > 0, while the corresponding rotating VA pair becomes vanishingly small. 
This result is of central importance for our main argument and will be analyzed in some detail in the 
remainder of this section. 

The best way to describe a vanishing VA pair, in the limit L — > 0, is to invoke yet another model 
configuration through the stereographic variable 



c- 



( + 



( = x + iy , C 



(32) 



where the constant d is taken to be real for simplicity. Now, configuration (|3"2"|) reaches a finite value 
fl = 1 at spatial infinity and the corresponding magnetization is uniform: 



m = (1,0,0) 



ICI 



(33) 



which is an appropriate boundary value (modulo a constant azimuthal rotation) for easy-plane anisotropy 
discussed here. Because of (|33[) the skyrmion number is expected to be an integer and is actually computed 
to be AC — 1 by a direct application of Eq. (fTU|) . Furthermore, the spin configuration derived from (f3"2")l 
is an exact solution of the LL equation, if we neglect anisotropy, with exchange energy 

' ^ (Vm ■ Vm) dxdy = An , (34) 



E = E P 



for any d, as discussed long time ago by Belavin and Polyakov [271 ]. 

We now return to the main line of argument by noting that configuration (|32p may also be thought 
of as a topologically nontrivial (AC = 1) model VA pair consisting of a (n, A) = (1,-1) vortex and 
a (k, A) = (—1,1) antivortex located at a distance d apart, in close analogy with the model VA pair 
constructed in Section [TTJ Although (|32p is not an exact solution in the presence of anisotropy, it provides 
a good model for the behavior of a vanishing VA pair in the limit L —* 0. Anisotropy sets a distance scale 
R ~ = 1 ( m rationalized units) beyond which a physically acceptable configuration must reach the 

uniform value ([33]) . In a sense, this property is shared by configuration (|32|) because it becomes uniform 
almost everywhere when d <ti 1 while it retains its topological structure as long as d remains finite. The 
strict limit d — ^ is not uniform. If taken naively, all topological structure appears to be lost and both 
energy E and skyrmion number Af appear to vanish. However, if integrals are performed before taking 
the d — * limit, E approaches 4-7T without encountering an energy barrier while AC = 1 for all d. Clearly 
the limit d — > creates a singular point which hides all topological structure. The main point is the claim 
that a similar situation arises in the calculated rotating VA pair in the limit L — > 0, as discussed further 
in Section fVl 

This section is completed with a comment concerning the manner in which the L — > limit is reached. 
Our numerical data as well as virial relation (|2T[) are consistent with a linear dispersion 



E 



(35) 



11 



in the limit L — > 0, which implies a finite value of the angular velocity u> — dE/dL = 1/2 in the limit of 
a vanishing VA pair, and thus an upper limit aj max ~ 1/2. 



V. VORTEX CORE SWITCHING 



Applied for long time intervals, our relaxation algorithm revealed some tendency for instability of a 
rotating VA pair, probably due to radiation effects analogous to those expected for a pair of rotating 
electric charges discussed in the Appendix. However, the basic features of rotating VA pairs studied 
in the preceding section persist over sufficiently long time intervals and are thus relevant for practical 
applications. In fact, in a realistic ferromagnet, some dissipation is always present and can be modeled 
by introducing Gilbert damping in the LL equation through the replacement 

dm dm ( dm\ . dfl ,. . 90 

^r^^r +a { mx ^r) or 1 m^ {l - a) m (36) 

in Eq. (TT]) or Eq. (jSJ), respectively, where a is a dissipation constant. 

The dynamics of a topologically nontrivial (Af = ±1) VA pair may thus be summarized as follows. 
The vortex and the antivortex rotate around each other, while the pair shrinks due to dissipation. The 
energy of the pair follows approximately the curve of the right panel of Fig. [7] as its size (and its angular 
momentum) decreases. At vanishing size a singular point of the type discussed in the preceding section 
would be created and the total energy would reach the finite value E — An (in rationalized units). 
However, the discreteness of the lattice actually interrupts the process when the size of the pair becomes 
comparable to the lattice spacing. The VA pair disappears (i.e., the skyrmion number changes abruptly 
from Af = ±1 to J\f = 0) and a burst of energy equal to An is released into the system, probably in the 
form of spin waves. In physical units, the amount of energy released is given by E — 8ntA where t is 
the film thickness and A the exchange constant. For typical values t = 10 nm and A = 1CP 11 J/m we 
obtain the estimate E ~ 2.5 x 10 -18 J which is apparently in rough agreement with numerical simulations 



The scenario described above explains how a topologically forbidden (AAf = 1) transition can take 
place in a real ferromagnet but does not by itself account for the experimentally observed vortex core 
switching. The complete scenario involves two distinct steps. First, application of a short burst of an 
alternating magnetic field creates a VA pair in the vicinity of a preexisting single vortex. Second, a 
three-body collision takes place during which a suitable pair of vortices is annihilated through a AAf = 1 
transition of the type described above, and the final product is a single vortex with polarity opposite 
to that of the original vortex [lj|. We also note that a system with three vortices (two vortices and an 
antivortex which form a cross-tie wall) was observed in a rectangular platelet in [30j j. Their spectrum of 
eigenmodes was studied experimantally. 

Here we do not address the question of how a VA pair is actually created. Rather we concentrate on 
step 2 of the process and explain in some detail how vortex core switching may occur in a three-body 
collision. Specifically, let us assume that a single (1, — 1)=C vortex is initially at rest at some specified 
point which is taken to be the origin of the coordinate system. Let us further assume that a topologically 
trivial (Af — 0) VA pair consisting of a (1, 1)=A vortex and a (—1, 1)=B antivortex is somehow created 
in the neighborhood of the original vortex. Once created the AB pair will undergo Kelvin motion of the 
type described in Section [TTT1 and eventually collide with the single vortex C. 

The process was simulated by a numerical solution of the corresponding initial- value problem in the LL 
equation. Figure [8] provides an illustration with three characteristic snapshots in the case of a relatively 
slow Kelvin pair initially moving along the y axis with velocity v = 0.1 for which the vortex and the 
antivortex are separated by a distance d ss 1/v — 10. As the pair approaches, the original C=(l,— 1) 
vortex teams up with the B=(— 1,1) partner of the AB pair to form a new, topologically nontrivial 
(Af = 1) VA pair in quasi-rotational motion. In fact, B rotates almost a full circle around C before 
rejoining its original partner A. The new AB pair is again a topologically trivial (Af — 0) VA pair in 
Kelvin motion that moves away from the target vortex, having suffered a total scattering angle that is 
greater than ir/2 from its original direction. The scattering is inelastic in the sense that the outgoing AB 
pair moves out with greater velocity (v = 0.15). And, most remarkably, the target vortex C moves away 
from the origin and comes to rest at a new location in the fourth quadrant of the xy plane. 

This rather unusual behavior is explained by the unusual nature of the conservation laws (11) and 
(12) which allow for a transmutation between position and impulse in the case of topologically nontrivial 
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FIG. 8: Three snapshots for the collision of a VA pair in Kelvin motion (the AB pair), initially located at (0, —15) 
and propagating with velocity v = 0.1, against a target vortex C initially located at the origin. During collision, 
antivortex B rotates around vortex C before rejoining its original partner A to form a new VA pair that scatters 
off at an angle in the third quadrant. The target vortex C is shifted to a new position in the fourth quadrant 
thanks to transmutation of VA pair momentum to vortex position. 



systems, such as the three-vortex system considered here (a more detailed discussion will be given in a 
future publication). 

The preceding numerical experiment was repeated for a Kelvin pair with relatively large velocity v — 0.5 
for which the vortex and the antivortex are tightly bound at a relative distance d = 2.6 [13] • The process 
is again illustrated by three characteristic snapshots in Figure^ While the initial stages of the process are 
similar to those encountered in the case of slow Kelvin motion (Figure [5]) a substantial departure occurs 
when the pair now approaches the target vortex. In particular, as soon as antivortex B=(— 1, 1) begins 
to rotate around the target vortex C=(l, —1) they collide and undergo a spectacular Aftf — 1 transition 
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FIG. 9: Same as Fig. [8] but for a larger initial velocity v = 0.5 of the AB pair. During collision, antivortex 
B begins to rotate around vortex C but the rotating BC pair is eventually annihilated leaving behind vortex A 
(with polarity opposite to that of the target vortex C) and a burst of spin waves that propagate away from the 
scattering region. 



(annihilation) leaving behind the A=(l, 1) vortex which may be thought of as the target vortex C=(l, —1) 
with polarity flipped from —1 to 1 (vortex core switching) and a burst of spin waves propagating away 
from the scattering region. 

A detailed numerical investigation of the three-vortex process for Kelvin waves with velocities in the 
allowed range < v < 1 suggests the existence of the three characteristic regions separated by two critical 
velocities v\ — 0.3 and v-i = 0.9 (such that < v\ < vq < «2 < 1, with vq = 0.78 being the critical 
velocity discussed in Section IIIip . For < v < v% , the Kelvin pair undergoes nearly elastic scattering of 
the type depicted in Figured For v\ < v < V2, the process leads to a topologically forbidden AM = 1 
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transition of the type illustrated in Figure [9l There is also some evidence that fast Kelvin waves with 
velocities in the narrow range v% < v < 1 undergo a nearly elastic scattering without inversion of the 
polarity of the target vortex. 



VI. CONCLUSION 



The VA pairs analyzed in this paper are special examples of solitary waves whose dynamics is closely 
related to their topological structure. For example, the VA pairs studied in Section Hill can undergo free 
translational motion because their topological charge vanishes (Af — 0). In contrast, a VA pair with 
nonvanishing Af performs rotational motion around a fixed guiding center and is thus spontaneously 
pinned within the ferromagnetic medium, as discussed in Section ITVl Such a peculiar dynamical behavior 
would have been surprising had it not occurred previously in the case of interacting electric charges in 
the presence of a magnetic field. 

It is then interesting to ascertain conditions under which a certain field theory would exhibit a similar 
link between topology and dynamics. A simple criterion was introduced in Ref. 22] and is briefly 
summarized as follows. We restrict attention to 2D Hamiltonian systems described in terms of A pairs 
of canonically conjugate variables (IF, $j) with i = 1, 2, ... A. Then one may define the local vorticity 

A 
i=l 

which is a simple generalization of the first step of Eq. (|9|), the remaining two steps being special to 
the specific example considered in the present paper (ferromagnets). Now, in general, the total vorticity 
r = J 7 dxdy is expected to vanish by a trivial partial integration using Eq. (|57)) . On the other hand, a 
nonzero F would signal a special topological structure of the field theory under consideration and may 
lead to peculiar dynamics. The theory analyzed in the present paper is an example (with A = 1) which 
yields a nonzero T that may be identified with the Pontryagin index Af — r/47r. An example with A = 2 
is provided by a 2D antiferromagnet where T = except when an external field is present which may 
lead to r 7^ and an interesting link between topology and dynamics [22| . 

Implicit in the preceding general argument is the fact that a topological charge Af is conserved. Also 
taking into account that Af is quantized, one would expect that a topological (Af ^ 0) soliton cannot be 
annihilated in a continuous manner. Nevertheless, a quasi-continuous process was described in Ref. [lj| 
and in the present paper according to which a rotating VA pair with Af =1 may be reduced to a singular 
point and thereby be eliminated by lattice discreteness without encountering an energy barrier. 

A mechanism for changing the topological number of a magnetic configuration makes it possible to 
obtain controlled switching between topologically distinct (and thus robust) magnetic states. This was 
achieved in the experiments of Refs. [13, [lj]. The dynamics underlying both experiments involves a 
three-vortex process [IH initiated by the production of a topologically trivial VA pair in the vicinity of 
a preexisting vortex by an alternating magnetic field. It has been predicted that the same phenomenon 



would occur if one uses a rotating external field 31| . The resulting three- vortex system carries nonzero 



topological charge and is thus by itself a rotating object spontaneously pinned in the magnet. Also due 
to dissipation, a quasi-continuous process takes place that changes the topological number by one unit, 
leaving behind a burst of energy in the form of spin waves and a single vortex with polarity opposite to 
that of the original vortex. 

Although our strictly 2D treatment provides a detailed scenario for the three- vortex process that leads 
to polarity switching, it does not account for the initial production of a topologically trivial VA pair. This 
is partly due to our approximation of a thin film with infinite extent. Implicit in this approximation is 
the assumption that the demagnetizing field amounts to a simple (additive) renormalization of easy-plane 
anisotropy. While this assumption appears to be firmly established for static magnetic states [19j | we do 
not know of a corresponding mathematical derivation for dynamical processes of the type discussed here. 

It is also important to visualize how the formation of a singular point discussed in the present strictly 
2D context appears within a realistic magnetic element of finite extent. Numerical simulations [T5 | show 
that, in an element of finite thickness, a singular point is first created at one of the surfaces of the element. 
The VA pair then vanishes by formation and subsequent annihilation of a singular point at successive 
levels away from the surface. At the stage when a singular point has been formed and annihilated, 
say, near the top surface, while the VA pair is still present in the bulk of the element, a Bloch Point 
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(BP) is created in the element. This is a somewhat simplified realization of the BP studied in Ref. [32[. 
Needless to say, the BP created near the top surface is eventually annihilated when the VA pair exits the 
system through the lower surface. It is important to emphasize the unusual fact that during creation and 
annihilation of the BP the system does not have to overcome an energy barrier, unlike the case discussed 
in [33J, essentially for the same reasons explained for the strictly 2D VA pairs studied in the main text. 
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APPENDIX A: ELECTRIC CHARGES IN A MAGNETIC FIELD 

Most of the distinct features of the dynamics of VA pairs occur also in the dynamics of electric charges 
in the presence of a uniform magnetic field B. Although we shall mainly be interested in 2D motion in a 
plane perpendicular to B, it is convenient to keep for the moment 3D notation and write the equations 
of motion for two interacting charges ei and e 2 : 

m^T =fi+ei(«i x B), m — ^ = F 2 + e 2 (v 2 x B) , (Al) 

where 

F l = -F 2 = - , ri " r2 , V'(\n - r 2 \) (A2) 
In - r 2 \ 

is the mutual force derived from a potential energy V = V(\ri — r 2 \) and V denotes derivative with 
respect to the argument \r\ — r 2 \. The conserved energy functional is then given by 

E = ±m(v 2 1 +v 2 2 ) + V(\r 1 -r 2 \), (A3) 

which does not depend explicitly on the magnetic field, while the conserved linear momentum (impulse) 
is now given by 

P = m(v! + v 2 ) - (ein + e 2 r 2 ) x B (A4) 

and differs from the usual mechanical definition by an important field dependent term. This (second) term 
actually indicates a rather profound influence of the magnetic field on the dynamics of electric charges. 
For instance, note that a shift of the origin of coordinates by a constant vector c, thus ri — > r± — c and 
r 2 — > t*2 — c, induces a nontrivial change on the impulse P of Eq. (|A4|) given by 



P —> P + {ei + e 2 )(c x B). (A5) 

This unusual behavior is analogous to that of the impulse defined by Eq. (fTTj) in the case of field config- 
urations with nonvanishing total topological vorticity T or skyrmion number TV. Here we may abstract 
from Eq. (|A5j) the analog of the topological vorticity in the present problem: 

T~(e 1 +e 2 )B. (A6) 

An electron-positron pair (r = 0) may undergo Kelvin motion, while two like charges (r 7^ 0) perform 
rotational motion around a fixed guiding center, as determined by the explicit solutions constructed and 
briefly analyzed in the remainder of this appendix. 

Consider first the case of an electron-positron pair (ei = —e 2 = e) for which a special 2D solution of 
Eqs. (|A1|) is given by 

d 

xi = - , yi = vt , z x = , 
x 2 = -~, y 2 = vt, z 2 = 0, (A7) 
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FIG. 10: Energy E and impulse P as functions of velocity v (left panel) and E vs P dispersion (right panel) for 
an electron-positron pair (ei = — ez) in Kelvin motion. 



where the electron and the positron are located at a constant relative distance d along the x axis and 
move in formation along the y axis with constant velocity 



V'(d) 
eB 



(A8) 



in close analogy with the Kelvin motion of the VA pair discussed in Section IIIII We further calculate 
energy E from Eq. (|A3|) and impulse P = (0, P, 0) from Eq. (|A4j) to find 



E = mi) 2 + V{d) = m 
P = 2mv + eBd = 2m 



V'(d) 
eB 

V'(d) 
eB 



+ V(d) 
eBd. 



(A9) 



Hence all relevant quantities are given in parametric form as functions of relative distance d once the 
potential energy V = V(d) has been specified. But there are some generic properties of this solution that 
are practically independent of the choice of V(d). We may take derivatives with respect to d of both sides 
of Eq. (jUJ to find E' = V'(l + 2mV"/e 2 B 2 ) and P 1 = eB(l + 2mV" / e 2 B 2 ). An immediate consequence 
of these relations is the group- velocity relation v = dE/dP. We also note that both E and P may acquire 
an extremum at a common value of d (or v) determined from 



1 



2mV"(d) 
e 2 B 2 



0. 



(A10) 



This is actually the reason for the appearance of a cusp in the E vs P dispersion analogous to that 
encountered in the Kelvin motion of VA pairs, which now appears to be a generic feature of a wide class 
of physical systems. 

For an explicit demonstration we make the special choice of potential energy 



V = 2ir In In - r 2 | 



(AH) 



in order to model the behavior of VA pairs at large relative distance 1(| 34 [ . For a graphical illustration 
we also make the special choice of constants m = 1 and eB ~ 2ir, as suggested by Eq. (|A6|) . to write 



1 

d' 



E = 



d 2 



27rlnd, P = 



2nd 



(A12) 



where we note that both E and P acquire a minimum at a distance d = do = 1/\/tt or velocity v = vq = 
y/n. The dependence of E and P on the velocity v as well as the E vs P dispersion are shown in Fig. fTOl 
and are found to be closely analogous to the results of Fig. [4] pertaining to Kelvin motion of VA pairs. In 
particular, for a widely separated pair (branch I) we find from Eq. (|A12[) that P ~ 2nd, vP ~ 2ir, v = 
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1/d, and E = 2ir ln(P/Po) with Po — 2tt, in close analogy with the asymptotic results of Eqs. (f22|) -(l24 |) . 
The appearance of a cusp and consequently of branch II in the spectrum is also notable. But the details 
of branch II arc different than those of Fig. [4] and Eq. (f26|) . There is now no upper limit in the velocity v. 
In fact, all v, E and P in Eq. (|AL2[) diverge in the limit of small d and E ~ P 2 /4, which coincides with 
the dispersion P 2 /2M of a free particle with mass equal to the total mass of the pair (M = 2m = 2). 

We have also carried out a stability analysis of the special Kelvin-like solution (|A7|I to find that the 
motion is marginally stable along branch I but becomes unstable along branch II. This conclusion is in 
agreement with a similar result obtained in Ref. [26] in the case of a vortex ring in a superfluid, as is 
further discussed in the concluding remarks of our Section IIIII 

As a last example we consider the case of 2D motion of two like charges e\ = e2 = e. Then a special 
solution of Eq. (|A1|) is given by 

x\ = —X2 = R cos cot , yx — — j/2 = R sincot , z\ = Z2 = (A13) 

which describes a pair rotating at constant radius R = d/2 and angular frequency lo = v/R where the 
velocity v is determined from the algebraic equation 



mv 2 



+ eBv - V'{d) = (A14) 



that expresses the exact balance of the centrifugal, the magnetic, and the mutual force. The conserved 
energy is still calculated from Eq. (|A3|) with v\ = v\ = v 2 , but the conservation of the impulse P of 
Eq. (|A4[) is simply equivalent to the statement that rotation takes place around a fixed guiding center. 
More relevant is now the conserved angular momentum (impulse) which is given by 

&\P> &2B 
L = m{xiyi ~ yi±i) + -^-(#1 + vl) + ™(^22/2 - 2/2±2) + —^~{ x \ + vl)i ( A15 ) 

where the overdot denotes time derivative. Again, the angular momentum differs from its standard 
mechanical expression by important field-dependent terms. 

Now, for the specific choice of potential energy given by Eq. (|A11[) . and constants m = 1, ei = €2 = e 
and eB = 2ir, the algebraic equation (| A14|) yields 

v = tt ^Jrz + I -r \ f R=^, (A16) 

while the angular velocity u>, the energy E, and the angular momentum L, read 

u = — , E = v 2 + 2tt ln(2i?) , L = 2Rv + 2n R 2 . (A17) 
R 

In view of Eq. (| A16|) all three quantities in (| A17|) are expressed in terms of a single parameter R = d/2. 
As a check of consistency one may verify the relation u> = dE/dL using Eqs. (|A16HA17|) . 
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The dependence of E and L on angular velocity u> as well as the E vs L dispersion are shown in Fig. [Til 
Again there exists a close analogy with the results of Section ITVl on rotating VA pairs. In particular, for 
large diameter d, Eqs. (|A17|) yield L ~ |d 2 , ujL ~ n, lu ~ 2/d 2 , and E 1 = 7rln(i/L ) with Lo = t/2, 
which should be compared with the asymptotic results for VA pairs given in Eqs. (|3U)) . On the other 
hand, some quantitative differences arise at small d, where the angular momentum vanishes as expected 
(L ~ y^nd) but the angular frequency diverges (w ~ 2^/n/d). Furthermore, the energy J5 does not reach 
a finite value at L — (as was the case for rotating VA pairs) but diverges logarithmically to minus 
infinity. 

Finally, an analysis of mechanical stability [35| shows that circular motion of two like charges in a 
magnetic field is marginally stable for all values of d, in contrast to the Kelvin motion discussed earlier in 
this section which becomes unstable at small d. However, unlike Kelvin motion which proceeds with no 
acceleration, a rotating pair is expected to radiate when full electrodynamics is turned on. Surely, this is 
also a source of instability and may indicate a similar instability for the rotating VA pairs discussed in 
Section [TV] 
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